Exotic p-wave superfluidity of single hyperfine state Fermi gases in optical lattices 
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We consider p-wave (triplet) pairing of single hyperfine state ultracold atomic gases trapped 
in quasi-two-dimensional optical lattices. We find that the critical temperatures in the lattice 
model is considerably higher and experimentally attainable around half-filling in contrast to the 
predictions of continuum model for p-wave superfluids. In tetragonal lattices, we show that the 
atomic compressibility and spin susceptibility have a peak at low temperatures exactly at the half- 
filling, but this peak splits into two in the orthorhombic lattices. These peaks reflect the p-wave 
structure of the order parameter for superfluidity and they disappear as the critical temperature 
is approached from below. We also calculate the superfluid density tensor, and show that for the 
orthorhombic case there is no off-diagonal component, however in the tetragonal case an off-diagonal 
component develops, and becomes a key signature of the exotic p-wave state. 
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Tunable optical lattices have been extensively used to 
study phase transitions in bosonic atomic gases |lj,|2|], 
cause they allow the controlled manipulation of the parti- 
cle density n, and of the ratio between the particle trans- 
fer matrix element t, and the interparticle interaction 
strength V [1 Q . This kind of control is not fully present 
in standard fermionic condensed matter systems, and 
has hindered the development of experiments that could 
probe systematically the effects of strong correlations as 
a function of n and t/V. However, fermionic atomic 
gases like 6 Li and 40 K have been succesfully trapped, 
and their normal state and superfluid properties are be- 
ginning to be studied |M S IE @- Since superfluid 
phases are more easily accessible in the experiments in- 
volving ultracold atomic gases, spin-polarized ultracold 
atomic systems are ideal candidates for the observation 
of novel triplet superfluid phases and for testing theoret- 
ical models that were proposed earlier. Thus, it is only 
natural to propose that optical lattices could be used to 
study the normal state and superfluid properties of ul- 
tracold fermionic systems as a function of n, t/V and 
lattice symmetry. This systems are of a broad interest 
not only for the atomic physics community but also for 
the nuclear, condensed matter and more generally for 
the many-body physics communities where superfluidity 
models have been investigated in various contexts. 

Presently there is only experimental evidence that K 
(Ref. HHGl) and 6 Li (Ref. [3 El El El El) 

can form 

weakly and tightly bound atom pairs, when the magnetic 
field is swept through an s-wave Feshbach resonance. 
However, the properties of p-wave spin-polarized ultra- 
cold fcrmions and their possible superfluid behaviour are 
beginning to be investigated [TR Il8j. When identical 
fermionic atoms are trapped in a single hyperfine state 
(SHS) the interaction between them is strongly influ- 
enced by the Pauli exclusion principle, which prohibits 
s-wave scattering of atoms in identical spin states. As 
a result, in SHS degenerate Fermi gases, two fermions 
can interact with each other at best via p-wave scatter- 



ing. Thus, one expects that the superfluid ground state 
of such SHS Fermi gases to be p-wave and spin triplet. 

In the p-wave channel, if the atom-atom interactions 
are effectively attractive then the onset for the formation 
of Cooper pairs in three dimensions occurs at a tempera- 
ture T c ss Ep exp[— 7r/(2fci?a sc ) 3 ] in the BCS limit, where 
Ep is the Fermi energy and a sc is the p-wave scattering 
length. Unfortunately, this temperature is too low to 
be observed experimentally. However, in the presence of 
the Feshbach resonances ^3, ^| , p-wave interactions can 
be enhanced, and the critical temperature for superfluid 
is expected to increase to experimentally accessible val- 
ues. On the other hand, we show that the spin triplet 
(p-wave) weak coupling limit in optical lattices (like in 
the singlet cases 8]) is sufficient to produce a superfluid 
critical temperature that is accessible experimentally. 

In this manuscript, we consider quasi- two-dimensional 
optical lattices with a periodic trapping potential of 
the form U(r) = J2iUo,i cos2 {hxi), with Uo, z ^> 
min{J7o,a:, C^o,j/}i which strongly suppresses tunneling 
along the z direction. This is a non-essential assumption, 
which just simplifies the calculations, but still describes 
an experimentally relevant situation. Here Xi = x, y, or z 
labels the spatial coordinates, ki = 2ir/\i is the wave- 
length, and t/o,i is the potential well depth along direction 
Xj, respectively. The parameters E/o.j are proportional to 
the laser intensity along each direction, and it is typically 
several times the one photon recoil energy Er such that 
tunneling is small and the tight-binding approximation 
can be used. 

Thus, in the presence of magnetic field h, we consider 
the following quasi-two-dimensional lattice Hamiltonian 
(already in momentum space) for an SHS Fermi gas 

h = J2 e(k)4 T «k T + \ E v ^ k ')Ok<, q , (!) 

k k,k',q 

where the pseudo-spin | labels the trapped hyper- 
fine state represented by the creation operator aj^, 

and b l,q = a k+q/2,T a -k+q/2,T- FurtnCmlOTC > £0) = 



2 



e(k) — Ji, describes the tight-binding dispersion e(k) = 
— t x cosk x a x —t v cosk y a y —t z cos k z a z , with fi = ^i+g/ish, 
and min-t^,^} > t z . 7(k,k') = uj t (k)Vw;(k / ) is 
the p-wave pairing interactions with matrix elements 
Vij = —2Vo^5ij and ii^(k) = (sinfc 2 ,a. E ,sinfc a a, / ). Here, 
Vo.i > is the effective interaction and a* is the corre- 
sponding lattice length along the i th direction. 

Using the functional integration formalism |l9l | (/? 
1 7 and units K = ks = 1), we obtain the saddle point 
effective action 



-_A — A 



where we use p 



(2) 



(p, ivi) with 



= (21 + l)?r//3 and 
define stationary vector field Aq = (Aq^A^). Here 

G- x (p)//3 = iu«7 -f (k)o- 3 + Ajw(k)a_ + wt(k)A cr + is 
the inverse Nambu propagator and <r± = (<7i ±<T2)/2 and 
<jj is the Pauli spin matrix. The condition SSq/SAq = 
leads to the order parameter equation 

A = MA (3) 

where M has matrix elements My = 
J2k Vo,i sin fcjaj sinfcjOj tanh(/3£'(k)/2)/£'(k). Here, 
E(k) = (£ 2 (k) + |A(k)| 2 )2 is the quasi-particle energy 
and the scalar order parameter A(k) = u;^(k)Ao is 
separable in temperature T and momentum k. 

Within the irreducible representations of the T>4h i^ih) 
group in the tetragonal (orthorhombic) lattices, |20j our 
exotic p-wave state corresponds to the 3 E u (n) represen- 
tation with a (Z-vector given by 



d(k) = /(k)(l,z,0), 



(4) 



where /(k) = AX + BY, and X and Y are smk x a x and 
sinkyCiy, respectively. Notice that, this state also breaks 
time reversal symmetry, as expected from a fully spin- 
polarized state. In the tetragonal lattice, the stable so- 
lution for our model corresponds to the case A = B ^ 0, 
and thus to the 3 E u (d) representation, where spin-orbit 
symmetry is preserved, but both spin and orbit symme- 
tries are independently broken. In the orthorhombic lat- 
tice, the stable solutions correspond to either A ^ 0, B = 
0or4 = 0,B^0, thus leading to the 3 E u (b) represen- 
tation. 

Our main interest is in tetragonal (square) lattices, 
however, we also want to investigate the effects of small 
anisotropies in optical lattice lengths. For defmiteness, 
we set a,y = a constant and investigate square and or- 
thorhombic lattices where a x = a and a x — a(l — 8) 
with <5 <C 1, respectively. In the case of square lat- 
tices, we choose the parameters t x = t y = t and Vq }X = 
Vo, y — Vo — 0.3i and change Vq jX and t x accordingly as 
we vary a x . Using exponentially decaying on-site Wan- 
nier functions and the WKB approximation, we obtain 
that the tunneling and interaction matrix elements along 



the rr-direction are proportional to exp (28WUo^x/Er). 
In the present calculation, we take t x = t exp (10(5) and 
Vq,x — Vo exp (10<5). As we increase (decrease) the ratio 5, 
both interactions and tunneling rate increase (decrease) 
in x direction. A different choice of on-site Wannier 
functions does not change our general conclusions, the 
only qualitative difference is that t x and Vo. x are differ- 
ent functions of 5. Depending on the lattice anisotropy, 
there are three distinct solutions: (1) Ao >x ^ 0, Aq )3/ = 
{A^0,B = 0); (2) Aq,, = 0, A , a ?0(A = 0,B? 0); 
and (3) A 0i a; = A , a / (4 = 5 / 0). In a square 
lattice both directions are degenerate (case 3), but even 
a small anisotropy in the lattice spacings or in the lattice 
potential lifts the degeneracy and throws the system into 
either case (1) or (2). 

The critical temperature, T c = max{T c ^ Xl T c ^ y }, is de- 
termined from the condition det M = 1 in Eq. (J3J, and 
can be written as 



°= n 



sin 2 kid 



£(k) 



-tanh^ k) 



2T r 



(•5) 



But, this equation has to be solved simultaneously with 
the saddle point number equation N = —dflo/djl where 
Oo = So(A(k))/{3 is the saddle point thermodynamic 
potential. This leads to N = X)k n (k) where 



n(k) = \ 



l^tanh^M 
E(k) 2 



(6) 



is the momentum distribution. These equations are 
approximately valid for all temperatures (T < T c ) in 
the limit of weak interactions (BCS limit). Plots of 
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FIG. 1: (Colour online) Plot of a) critical temperature T r = 
T c /Eq, and b) chemical potential /i r = p,/Eo versus filling JV C 
for the square case a x = a (red); and the orthorhombic case 
a x = 0.99a (blue). 

T r = T c /Eq and fx r = fi/Eo as a function of number 
of atoms per unit cell (N c ) are shown in Fig. 1 for the 
cases of square and orthorhombic lattices. Here, E = 2t 
is the half-filling Fermi energy for the square lattice. No- 
tice that T r is maximal at half-filling, and that it has a 



3 



value 0.01, which is much higher than the theoretically 
predicted T c from the continuum model, and compara- 
ble to experimentally attainable T/Tp sa 0.01. This im- 
plies that the superfluid regime of spin-polarized fermion 
gases may be observed experimentally in a lattice, even 
in the limit of weak interactions. The observability of 
a superfluid transition in spin-polarized fermion systems 
is clearly enhanced when the system is driven through a 
Feshbach resonance (in a lattice or in the continuum), as 
T c is expected to increase further in this case, however 
our calculations indicate that the weak interaction (BCS) 
limit may be sufficient in the lattice case. 

Notice that the order parameter symmetry dramati- 
cally effects T c as can be seen by rewriting Eq. J^J as 

„ rt*+t tanh 

1 = J] V ^ / d£ of : D pA^ (?) 

i=x,y J -^~ t 

where we define an effective density of states (EDOS) 
D Pt i{e) — ^2^S[e — e(k)](v2sin /ciai) 2 which is plotted 
in Fig. 2 for a x = a (black), a x = 0.99a (red), and 
a x = 0.95a (blue). Here, \^2smkiai are symmetry fac- 
tors related to the order parameter. For instance in the 
s-wave case, this symmetry factor is 1 and EDOS be- 
comes D{e) — J2k^i £ ~~ e (k)L which is the density of 
states (DOS) of normal fermions. For the p-wave symmc- 
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FIG. 2: (Colour online) Plot of a) normal state DOS D(e) 
and b) p-wave EDOS D P:X (e) versus energy e/Eo in units of 
t for a x — a (black) and a x = 0.99a (red), and a x = 0.95a 
(blue) in two dimensions. 

try discussed, D p ^(e)Vo^ plays the role of a dimcnsionlcss 
coupling parameter which controls the critical tempera- 
ture. Notice that, EDOS and T c are maximal at half- 
filling in square lattices (a x = a). Additionally, EDOS 
decreases and finally vanishes at the band edges where a 
small ratio of T c /Tp was predicted from continuum mod- 
els. However, this is not the case around half-filling and 
we expect weak interactions to be sufficient in observing 
superfluidity. In the orthorhombic lattices, EDOS and T c 
are considerably increased and are maximal around half- 
filling with a small anisotropy (a x = 0.99a) in the lattice 
spacings. Notice, however, that further anisotropy in the 



lattice spacings (a x = 0.95a) leads to a decrease in EDOS 
and T c around half-filling with respect to the case with 
a x — 0.99a. 

Next we discuss the atomic compressibility and the 
spin susceptibility of the system. The isothermal atomic 
compressibility is given by k = — {\/N 2 )d 2 £l/dji 2 — 
{\/N 2 )dN c /djl where 

dfL ~ 2£3(k) 2 +l. r(k) £2 (k) (8) 

and Y(k) = (/3/4)sech 2 [/3E(k)/2]. Plots of n r = k/k q , 
where normalization kq is evaluated at T = and half- 
filling, is shown in Fig. 3 for the square and orthorhombic 
lattices. In the square lattice case, k has a peak at half- 
filling and low temperatures, and a hump at T c . The 
same qualitative behavior is obtained in the case of or- 
thorhombic lattices, with the additional feature that the 
central peak (hump) splits into two due to degeneracy 
lifting. 
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FIG. 3: (Colour online) Plot of compressibility k t = k/kq 
versus filling N c at T — 0.004T c maa; for a) a x = a (red) and 
b) a x = 0.99a (blue). Notice that the compressibility peaks 
disappears at T = T c and turns into humps (black) in both 
cases. 

The peaks at T « 0, can be understood by noting 
that N 2 k can be written as ^ k g(k) / E(k) for T = 0, 
where g(k) — 2n(k)[l — «(k)]. Therefore, the peaks are 
due to non- vanishing <?(k) in regions of k-space where 
E(h) vanishes and n(k.) is rapidly changing. In tetrago- 
nal lattices, the integrand <?(k) / E(k) has 4 k-space points 
(0,±7r), and (±7r,0) in the first Brillouin zone (1BZ) 
where it diverges only when the chemical potential fi = 
(N c — 0.5). Similarly in orthorhombic lattices, the in- 
tegrand diverges at 2 k-space points (±7r, 0) in the 1BZ 
when fj, = t x -t = 0.021 (N c = 0.548). Furthermore, 
the integrand diverges at 2 k-space points (0, ±7r) when 
fi = -t x + 1 = -0.021 (N c = 0.452). For every other fi 
the integrand is well-behaved resulting in a smooth k in 
both cases. At T = T c , the humps are not related to the 
order parameter, but are due to the peaks appearing in 
DOS (Fig. 2a). Notice that, while DOS have only one 
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divergence at half-filling in the tetragonal case, it has two 
peaks in the case of orthorhombic lattices which leads to 
two humps. 

Furthermore, for a magnetic field h = hfj applied along 
an arbitrary 77-direction, the spin susceptibility tensor 
component is \vv = —d 2 Q/dh 2 = g 2 p: 2 B dN c j 'djl. There- 
fore, in spin-polarized systems, Xnv is directly related to 
the atomic compressibility and is given by [N 2 / (<7/j,s) 2 ]/t. 

Next we discuss the superfiuid density of the system. 
Generally speaking, there are two components to the re- 
duction of the superfiuid density at non-zero tempera- 
tures, one coming from fermionic (quasi-particle excita- 
tions) and the other bosonic (collective modes) degrees 
of freedom. The fermionic component comes from phase 
twists of the order parameter and the temperature de- 
pendence of its components is given by 



Pij 



(9) 



where n(k) is the momentum distribution, and di de- 
notes the partial derivative with respect to Here, we 
will not discuss the bosonic contribution, except to say 
that at low temperatures the dominant terms come from 
Goldstone modes associated with the phase of the order 
parameter, which are underdamped in our case due to 
sub-critical Landau damping. Furthermore, in the case 
of tetragonal symmetry, Goldstone modes do not con- 
tribute to the off-diagonal component of the superfiuid 
density, which is the main focus of the analysis that fol- 
lows. 
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FIG. 4: (Colour online) Plot of diagonal {p%* = p v r v \ upper 
curves) and off-diagonal (p% v = p^, x ; lower curves) superfiuid 
density tensor components p 1 ^ — pijj p™y X versus temperature 
T/T c at fillings iV c = 0.5 (black), 0.45 (red) and 0.4 (blue) 
when a x = a. 



In Fig. 4, we plot p l J = Pijj p™y X as a function of tem- 
perature for three fillings; N c = 0.5 (black), iV c = 0.45 
(red), and iV c = 0.4 (blue). Normalization p™y X is the 
maximum value of the off-diagonal component and is 
evaluated at half-filling and T = 0.02T C in a square lat- 
tice. It is important to emphasize that square lattices 



have identical diagonal elements p xx = p vv due to the 
tetragonal symmetry, but have a nonzero p xy component 
(Fig. 4) as a result of the absence of reflection symmetry 
in the yz-plane [x — > —x) or the xz-plane (y — > —y) for 
the d- vector defined in Eq. Q. Notice that, reflection 
symmetry is restored and p xy vanishes identically in the 
orthorhombic case for any temperature T. We show that 
while p xy is zero at T = for any filling, it increases 
until all tensor components have the same value at some 
higher temperature. 

In summary, we considered p-wave pairing of single- 
hyperfine-state Fermi gases in quasi-two-dimensional op- 
tical lattices. We found that the critical temperatures 
in tetragonal and orthorhombic optical lattices are con- 
siderably higher than the continuum model predictions, 
and therefore, experimentally attainable. At low tem- 
peratures, we found a peak in the atomic compressibility 
(and similarly in the spin susceptibility) exactly at half- 
filling for the tetragonal lattice. This peak splits into two 
smaller peaks in the orthorhombic case. These peaks 
reflect the p-wave structure of the order parameter at 
low temperatures, and they decrease in size as the criti- 
cal temperature is approached from below. Furthermore, 
in the orthorhombic lattices, the off-diagonal component 
p xy of the superfiuid density tensor vanishes identically, 
while the diagonal components p xx and p yy are different. 
However, in the square lattices, we showed that p xy 7^ 0, 
while p xx = pyy. The presence of non-zero p xy is a key 
signature of our exotic p-wave triplet state. 

Acknowledgement: We would like to thank NSF 
(DMR-0304380) for financial support. 



[1 
[2 
[3 
[4 
[5 
[6 

[? 

[8 
[9 

[10 

[11 
[12 
[13 
[14 
[15 
[16 
[17 
[18 
[19 

[20] 



M. 
M. 
C. 
K. 
M. 
G. 
D. 
W 
M. 



Greiner et al., Phys. Rev. Lett. 87, 160405 (2001) 
Greiner et al., Nature 415, 39 (2002) 
A. Regal and D. S. Jin, Nature 424, 47 (2003) 
M. O'Hara et al., Science 298, 2179 (2002) 
Kohl et al. , | |cond-mat /0410389 
Modugno et al, Phys. Rev. A 68, 011601(R) (2003) 
Jaksch et al, Phys. Rev. Lett. 81, 3108 (1998) 
Hofstetter et al, Phys. Rev. Lett. 89, 220407 (2002) 
Rodriguez and P. Torma, Phys. Rev. A 69, 041602(R) 
(2004) 

C. A. Regal et al, Phys. Rev. Lett. 92, 040403 (2004) 
M. Greiner et al, Nature (London) 426, 537 (2003) 
K. E. Strecker et al, Phys. Rev. Lett. 91, 080406 (2003) 
M.W. Zwierlein et al, Phys. Rev. Lett. 91, 250401 (2003) 
T. Bourdel et al , |cond-m at /040309T1 
M. Bartenstein et al, Phys. Rev. Lett. 92, 120401 (2004) 
J. Kinast et al, Phys. Rev. Lett. 92, 150402 (2004) 
J. L. Bohn, Phys. Rev. A 61, 053409 (2000) 
C. A. Regal et al, Phys. Rev. Lett. 90, 053201 (2003) 
V. N. Popov, Functional integrals and collective excita- 
tions, Cambridge University Press, New York (1987) 
J. F. Annett, Adv. Phys. 39, 83 (1990) 



